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Abstract 

We present an ah initio density functional theory (DFT)-based study of h-BN 
domain size effect on band gap of mono-layer li-BNC lieterostructure mod- 
eled as (B3N3)^(C6)i-z. The atomic structures, electronic band structures, 
density of states and electron localization functions of h-BNC are examined 
as h-BN concentration ranged from to 100%. We report that the electronic 
band gap energy of h-BNC can be continuously tuned in full range between 
that of two phases, graphene and h-BN, as a function of h-BN concentration. 
The origin of the tunable band gap in these heterosturcutes are due to the 
change in the electron localization with h-BN concentration. 

Keywords: Mono-layer Hexagonal BNC Heterostructures, Density 
functional theory. Band structure. Tunable band gap 



1. INTRODUCTION 

Many approaches have been proposed to impro ve the semicond ucting 



properties of g raphene including the use of strains ( iGui et al.l . l2008l ). pat 



terned defects ( lAppelhans et al.l . 120101 ) , dimensionality reductio n as nanorib- 



bons external supper-lattice potentials 
application of external electri c fields ( 



Ohta et al. 



cal species such as hydrogen ( lElias et al.l . 120091 : ISingh and Yakobsonl . l2009l ). 



Tiwari and Stroudl . l2009l ). doping and (or) 



2006). addition of chemi- 
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substrate surface chemistry (IShemella and Navakl. |2009| ) . and introduction 
of antidote lattice structures ( jPedersen et al.l . l2008l ) . Among all these meth- 
ods, one particularly promising technique has been reported to be the doping 
of B and N atoms into the graphene lattice. Experimental and theoretical 
studies show that both p-type and n-type semi-conducting graphene may be 
generated by subst ituting the C atoms with B and N atoms, respectively 
(IWang et al.l . |2009[ ). In such methods the dopant B and N atoms serve to 



modify the electronic structure of the graphene and increase the band gap. 

Hexagonal boron nitride (h-BN) and graphene have similar 2D lattice 
structures but with very different physical properties. Graphene has a zero 
band gap while h- BN is an insulator with a wide band gap (5.5 eV in bulk 
fISong et al.lJ20HjVL_ Interesting structure can be made by mixing these two 



structures (jSuenaga et al.l . 119971 ). It has been suggested that inserting a 



graphene layer between layers of h exagonal BN (h - BN) can exhibit tun- 



able band gap as large as 0.23 eV (ISlawihska et al.l . l2010l ). A small band 



gap can be also opened by placing a graphene layer onto a h-BN substrate 



(iGiovannetti et al.l . l2007l ) 



Recently, a promising method has been reported in which atomic mono- 
layers have been generated consisting of h-BN phases in gr aphene (h-BNC ) 
using a thermal catalytic chemical vapor deposition method ( ICi et al.l . |2010| ). 
This hybrid mono-layer has been shown to have isotropic physical properties 
which can be tailored by controlling the kinetic factors affecting the h-BN 
domain size within graphene. This is different from B-doped or N-doped 
graphene, where the integrity of the h-BN structure is missing. Although 
there are a few stud ies about domain- wise phase segregation in such B-C-N 



graphene-like layers (iPanchokarla et al.l . 120091 : IXu et al.l . l2010l : iManna and Pati 



201ll : lMartins and Chachaml . l201ll ). a systematic study of the electronic band 



gap with respect to the h-BN domain size is missing. 

The method of Adding h-BN phase is much more efficient in opening band 
gap in graphene than all other methods. There are two most important 
advantages in this method: (1) wide range of tunable band gap of up to 
4.7 eV can be obtained, which is one order of magnitude larger than other 
methods; (2) integrity of the individual phases is kept which make it much 
easier for fabrication. In this paper, we examine the h-BN domain size effect 
on band gap of h-BNC monolayer using ab initio density functional theory. 
The electronic band gap with respect to the h-BN domain concentrations is 
studied and the origin of the tunable band gap in these heterosturcutes is 
also discussed. 
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2. MODEL 



The h-BN domain in h-BNC hybrid structures is modeled by the con- 
centrations of h-BN while maintaining the hexagonal structure of BN locally 
within the system. The h-BN domain size effect then can be represented 
by the h-BN concentration x in the model as (B3N3)a;(C6)i-x where (B3N3) 
and (Ce) denote the nanodomain structure of h-BN monolayer and graphene, 
respectively. The proposed domain size effect model ('B3N3).-,;(' Cfi)i_.r is base d 



on th e results of previous stu dies of B.t;C,,Nz layered s t ructu r es (iMazzoni et al 



2006 ). layers and nanotubes dMartins and Chachanj. I2OIII : iBhowrnick et al 



201 ll ). quantum do ts and nanorods ( Singh and Yakobsonl . |2009| ). and mono- 



201l[ ). It is a general believe that there 



layer nano hybrids (iManna and Patil . 

are domain-wise segregation phases of h-BN and graphene Avith lower energy, 

larger band gap, and better t hermal stability ( ICi et al.l . l2010l : lMartins and Chachanj . 
201lUManna and Patil . lioilh . As a result, the six- atom hexagonal structures 
(both B3N3 and Cq) are the basic blocks in these heterostructures. This 
model of (B3N3)2.(C6)i-x captured the main feature of these heterostruc- 
tures. It was found that both stoichiometry and geometry vary the band 
gaps of B-C-N nanotube and l ayered materials and the maximum band gap 
is ach ieved at B /N ratio of 1 flMazzoni et al.l . 120061 : iMartins and Chachanj . 
20111 ). while in B-C-N monolayer, domain size is a dominant factor com- 



paring with domain shape in tuning band gaps (IManna and Patil . 1201 ll ) . By 
varying h-BN concentration x, the domain size effect on mechanical prop- 
erties of mono-layer hexagonal BNC heterostructures can be appropriately 
presented in our (B3N3)2.(C6)i-x model. 

We need to emphasis that although we used stoichiometry of h-BN here, 
this domain model is different from point model of B^C^N^ mo del at x = z 
wher e hexagonal (B3N3) structure is not considered, for example ( Mazzoni et al.l . 



20061 ). In other words, our model specified not only the stoichiometry of B/N 
ration of 1, but also the hexagonal (B3N3) and Cq structures, to represent 
the two separated phases in heterogeneous h-BNC structures. 

We examined the variation of elastic band gaps of the /i-BNC monolayer 
as function of /i-BN concentration. Five /i-BNC configurations, in order 
of /i-BN concentration, 0%, 25%, 50%, 75% and 100% have been studied, 
where 0% and 100% corresponds to the pure graphene and /i-BN, respectively. 
The three other concentrations were selected based on their simplicity and 
representativeness. The atomic structures of these five configurations (Fig. [T]) 
were determined by the ab initio density functional theory. 
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Due to the segregation of two phases h-BN (B3N3) and graphene (Cr) , 
the minority will form "island", as observed in experiments ( ICi et al.l . |2010| ). 
For a certain stoichiometry, as presented by h-BN concentration x, the larger 
system size used in the simulation box, the larger "island" will be formed 
to achieve lower energy and become more stable. In this study, we used a 
small system of 24 atoms due to its simplicity while containing its distinguish 
characteristics of modeling separated phases in h-BNC heterostructures. The 
system size effect of this model was examined in a larger system with 150 
atoms. 



3. COMPUTATIONAL DETAILS 

DFT calculations were carried out w ith the Vienna Ab-initio Simula- 
tion Package (VASP) f lKresse et al.lll993|l^ which is based on the Kohn-Sham 



Density Functional Theory (KS-DFT) (IHohenberg and Kohru . Il964j ) with the 



generalized gradient approximations as parameterized by Perdew, Bu r ke an d 



Ernzerhof (PBE) for exchange-correlation functions (iPerdew et al.l . Il996l ). 
The electrons explicitly included in the calculations are the {2s^2p^) electrons 
of carbon, the {2s'^2p^) electrons of boron and {2s'^2p^) electrons of nitro- 
gen. The core electrons (Is^) of carbon, boron and nitrogen are repl aced by 



the p rojector augmented wave (PAW) and pseudopotential approach (IBlochll . 



19941 ). A plane- wave cutoff of 520 eV is used in all the calculations. 

The criterion to stop the relaxation of the electronic degrees of freedom 
is set by total energy change to be smaller than 0.000001 eV. The optimized 
atomic geometry was achieved through minimizing Hellmann-Feynman forces 
acting on each atom until the maximum forces on the ions were smaller than 
0.001 eV/A. 

The atomic structures of the five configurations were obtained by fully 
relaxing a 24- atom-unit cell where all atoms were placed in one plane. The 
irreducible Brillouin Zone was sampled with a Gamma-centered 19 x 19 x 1 
A;-mesh and initial charge densities were taken as a superposition of atomic 
charge densities. There was a 14 A thick vacuum region to reduce the inter- 
layer interaction to model the single layer system. 

4. RESULTS AND DISCUSSION 

We examined the variation of band gap of h-BNC mono-layers as function 
of h-BN concentration. Five h-BNC configurations have been studied with 
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(a) (b) (c) . 




Figure 1: Atomic structures of five configurations in order of h-BN concentration: (a) 0%; 
(b) 25%; (c) 50%; (d) 75%; (e) 100%. (f) shows the Brillouin zone and high symmetrical 
k points. 



0%, 25%, 50%, 75% and 100% concentration of h-BN, with 0% and 100% 
corresponding to pure graphene and h-BN, respectively. The three other 
concentrations were selected based on their simplicity and representativeness. 
The atomic structures of these five configurations (Fig. [1]) were determined 
by the ah initio density functional theory. 

The optimized atomistic structures of the five configurations shown in 
panels (a)-(e) in Fig. [1] were determined using DFT calculations. The lattice 
constants are listed in Table [H Due to the intrinsic difference between 
pure h-BN and graphene, the lattice constants of the h-BNC mixtures are 
obtained by averaging the lattice vectors of the super-cells. We found that 
the lattice constant increases with h-BN concentration x. Our results are in 



good agreement wi th experiments in h-BN (2.51 A) ( iLiu et al. . 2003 ) and 
graphene (2.46 A) (jBaskin and Meyer . 1955 ). 

The electronic band structures of the optimized atomistic structures of 
the five configurations are plotted in Fig. [2]a-e, in order of x. The F, M and 
K are the high geometrically symmetrical k points in the Brillouin zone, as 
noted in panel (f) in Fig. [H The dashed line in the band structure plots 
notes the Fermi energy of graphene, which is the intrinsic Fermi energy for 
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Figure 2: Electronic band structures of five configurations in order of h-BN concentration: 
(a) 0.0; (b) 0.25; (c) 0.5; (d) 0.75; (e) 1.0. (f) shows symmetry of band structure in 
configuration (c). 
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Table 1: Lattice constant a, bottom edge of conduction bands Ec, top edge of valence 
bands band gap energy Eg of the five configurations from DFT calculations. 



X 


0.0 


0.25 


0.5 


0.75 


1.0 


a (A) 


2.468 


2.484 


2.496 


2.506 


2.512 


Ec (eV) 





0.57 


0.97 


1.72 


2.90 


(eV) 





-0.44 


-0.70 


-1.15 


-1.78 


E, (eV) 





1.01 


1.67 


2.86 


4.68 



h-BNC structures. The band energies are shifted so that the intrinsic Fermi 
surface energy is ahgned up. The bottom edge of conduction bands Ec, top 
edge of valence bands E^j, and band gap energy Eg of the five configurations 
from DFT calculations are also obtained and listed in Table [H We notice 
that the band gap increases with x, with the curvature of the band at F 
decreasing with x indicating increment of effective mass. The symmetry of 
the real space is broken in h-BNC since the elements are heterogeneous. We 
tested the symmetry of the electronic band structure at the point of M2 
and K2, which are center symmetrical points of M and K respectively. We 
found that the symmetry of electronic band structures at these points are 
still preserved, as shown in Fig. [2]F, where the band structure of configuration 
(c) is presented as an example. 

The density of states (DOS) of the optimized atomistic structures of the 
five configurations are plotted in Fig. [3] in order of increasing h-BN concen- 
tration. The lines were shifted by 40 (states/eV) to avoid overlapping. The 
energies were shifted so that the Fermi surface energy could be set to zero. 
This DOS plot clearly shows the systematic increment of the band gap with 
respect to h-BN concentration. 

The band gap energy (Eg) of the optimized atomistic structures of the five 
configurations from DFT calculations are plotted in Fig. H] as a function of 
X. Our calculation (DFT-GGA) of the Eg of graphene-lik e h-BN (4.67 eV) is 



in go od agreement with DFT-LDA prediction of 4.7 eV ( iGiovannetti et al. 



20071 ). What is most interesting is that Eg increases with x in a nearly 
quadratic manner. The relationship between Eg and h-BN concentrations 
may be approximated by Eg{x) = 2.16x^ -|- 2.42x as plotted in dashed line in 
Fig.H 

The system size effect to the band gap is examined with a 150-atom 
system of which 54 atoms in h-BN phase form an "island" surrounded by 96 
atoms in graphene phase as shown in Fig. |5l The band gap energy of this 
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Figure 3: Density of states of five configurations in order of h-BN concentration. 
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Figure 4: Band gap energy Eg as a function of h-BN concentrations x, fitted to a quadratic 
function. 
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54/150 h-BNC system is 1.001 eV as the diamond dot in Fig. IH agreeing 
well with the prediction from 24-atom model. The full density of states is in 
plotted in Fig. [5] 



Density of States 




15 -10 -5 
Energy (eV) 



Figure 5: The configuration (left) and Density of States (riglit) of the 54/150 (54 h-BN 
atoms in 150 h-BNC atoms) system. 



A continuous and wide range of tunable band gap of h-BNC materials is 
possible ranging from semi-metal (graphene) to semi-conductor {Eg up to 4 
eV) and insulator [Eg > 4 eV). Interestingly, the band gap energies of h-BNC 
at X = 0.25,0.5,0.75 (corresponding to BNCe, BNC2, B3N3C2) are similar 
to those of Sihcoi i (1.11 eV), Selenium (1.74 eV ) and Sihcon carbide (2.86 
eV), respectively ( Streetman and Banerjee . 1999[ ). which indicates potential 
substitutions for those semiconductors. 

The mechanism of this band gap energy increment may be roughly ex- 
plained with bond stretching caused by introducing h-BN in graphene. Since 
the lattice constant a increases with the h-BN concentration, the average 
bond length increases, with a reduction in the overlap of the wave functions, 
which causes the band gap enlargement. 

We also studied the projected density of states (PDOS) of the atoms in 
the system. Each atom in the system has the same band gap in PDOS as 
that of total system. The the atoms on the interface of h-BN and graphene 
were plotted in Fig. [61 Only Pz electrons, which form tt bonds in z direction, 
contribute to the band gap. The edge of the band gap becomes sharper 
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when the concentration of h-BN is higher, indicating the higher locahzation 
of electrons around the Fermi surface. 
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V % # % # 
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Figure 7: Electron localization function in the five configurations. 



In particular, since the total number of electrons of a h-BNC system is 
the same as in a corresponding graphene structure, the mechanism of this 
band gap tuning may be ascribed to electron redistribution, bonding and 
localization due to many-body interactions in the system. The "electron lo- 
calization function " (ELF) measures the extent of spatial localization of the 
reference electron ( Becke and Edgecombe . 19901 ). with values between and 
1. Values of 0, 0.5 and 1.0 correspond to vacuum, uniform electron gas and 
perfect localization, respectively. A higher ELF value at a point indicates 
that the electrons are more localized there than in a uniform electron gas of 
the same density. ELF plots of the five configurations from DFT calculations 
are shown in Fig. [7] a-e. The plots show that electrons are more localized 
around nitrogen atoms, which is consistent with the fact that nitrogen atoms 
have more valence electrons than carbon and boron atoms. The electrons on 
the B-N and B-C bonds are more localized than the C-C and N-C bonds. 
Since the increment of h-BN concentration enhances overall electron local- 
ization, the band gap energy enlarges as a consequence. 



5. CONCLUSIONS 



In summary, we used ah initio density functional theory to investigate the 
effect of the h-BN domain size on the band gap of h-BNC hybrid structures 



12 



with (B3N3)j.(C6)i_j; model. This simple model of h-BNC heterostructure 
with phase segregation not only considers the stoichiometry, but also the 
shape of the nano-structure: which should be hexagonal ring of (B3N3) or 
(Ce). The atomic structures, electronic band structures, density of states, 
projected density of states, electron localization functions of five h-BNC hy- 
brid structures were explicitly examined. We found that the band gap energy 
of h-BNC can be continuously tuned in full range between that of two phases, 
graphene and h-BN, as a function of h-BN concentration. The origin of the 
tunable band gap in these heterosturcutes are due to the change in electron 
localization by changing the h-BN concentration. This result may provide a 
new approach in practical engineering applications of these heterostructures. 
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